Shell model method for Gamow- Teller Transitions in heavy, deformed nuclei 
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A method for calculation of Gamow- Teller transition rates is developed by using the concept of 
the Projected Shell Model (PSM). The shell model basis is constructed by superimposing angular- 
momentum-projected multi-quasiparticle configurations, and nuclear wave functions are obtained 
by digonalizing the two-body interactions in these projected states. Calculation of transition matrix 
elements in the PSM framework is discussed in detail, and the effects caused by the Gamow- Teller 
residual forces and by configuration-mixing are studied. With this method, it may become possible to 
perform a state-by-state calculation for /3-decay and electron-capture rates in heavy, deformed nuclei 
at finite temperatures. Our first example indicates that, while experimentally known Gamow- Teller 
transition rates from the ground state of the parent nucleus are reproduced, stronger transitions 
from some low-lying excited states are predicted to occur, which may considerably enhance the total 
decay rates once these nuclei are exposed to hot stellar environments. 

PACS numbers: 21.60.Cs, 23.40.Hc, 23.40.-s, 23.40.Bw 
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I. INTRODUCTION 



The knowledge on weak interaction processes is one 
of the most important ingredients for resolving astro- 
physical problems. The first systematical work on stellar 
weak-interaction rates was performed by Fuller, Fowler, 
and Newman 0, [3.0. [| who recognized the decisive role 
played by the Gamow- Teller (GT) transitions. Due to 
their pioneering work, the study of stellar weak interac- 
tion rates requested by astrophysics becomes essentially 
a nuclear structure problem, in which the actual decay 
rates are determined by the microscopic inside of nuclear 
many-body systems. It has been suggested that the nu- 
clear shell model, i.e. a full diagonalization of an effective 
Hamiltonian in a chosen model space, is the most prefer- 
able method for GT transition calculations. This was 
noticed early by Aufderheide et al. and has recently 
been emphasized by Langanke and Martinez-Pinedo [(|. 

For a theoretical model employed in GT transition cal- 
culations, it is generally required that the model can re- 
produce a wide range of structure properties of relevant 
nuclei. It has been shown that the state-of-the-art shell- 
model diagonalization method is indeed capable of per- 
forming such calculations. For example, Wildenthal and 
Brown 0, Hj obtained nuclear wave functions in the full 
sd-shell model space, which were successfully applied to 
calculation of GT rates in the sd shell nuclei 9] . By using 
the method developed by the Strasbourg-Madrid group 
Langanke and Martinez-Pinedo 

El 

made the shell- 
model GT rates available also for the pf shell nuclei. 
Still, these sophisticated calculations are tractable only 
for nuclei up to the mass-60 region, and cannot be ap- 
plied to heavier nuclei which play important roles in the 
nuclear processes in massive stars. 



In the long history of the nuclear shell-model develop- 
ment, tremendous effort has been devoted to extending 
the shell-model capacity from its traditional territory to 
heavier shells. Despite the great progress made in recent 
years, it seems impossible to treat an arbitrarily large 
nuclear system in a spherical shell model framework due 
to the unavoidable problem of dimension explosion. One 
is thus compelled to seek judicious schemes to deal with 
large nuclear systems. The central issue has been the 
shell-model truncation. There are many different ways 
of truncating a shell-model space. While in principle, 
it does not matter how to prepare a model basis, it is 
crucial in practice to use the most efficient one. In this 
regard, we recognize the fact that except for a few lying 
in the vicinity of shell closures, most nuclei in the nuclear 
chart are deformed. This naturally suggests for a shell 
model calculation to use a deformed basis to incorporate 
the physics in large systems. That is the philosophy that 
the Projected Shell Model (PSM) Eland the several 
important generalizations 0, IT3. Il5l Il6t are based 
on. 

The present article reports on the new development for 
calculation of GT transition rates in the PSM framework. 
Before a detailed description of the work, we mention a 
few attractive features in our approach, which may be 
relevant for future astrophysical applications. 

• The PSM utilizes single particle bases generated by 
deformed mean-field models yet carries out a shell- 
model diagonalization like the conventional shell 
model. Conceptually, the PSM bridges two im- 
portant nuclear structure methods: the deformed 
mean-field approach and the conventional shell 
model, and takes the advantages of both. On the 
one hand, as a shell model, the PSM can be applied 
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to any heavy, deformed nuclei without a size lim- 
itation. On the other hand, unlike the mean-field 
models or models with an average nature, the PSM 
wave functions contain correlations beyond mean- 
field and the states are written in the laboratory 
frame having definite quantum numbers such as 
angular-momentum and parity. These are needed 
properties when the wave functions are employed 
in transition calculations. 

• Because of the way the PSM constructs its basis, 
the dimension of the model space is small (usually 
in the range of 10 2 — 10 4 ). With this size of basis, a 
state-by-state evaluation of GT transition rates is 
computationally feasible. This feature is important 
because in stellar environments with finite temper- 
atures, the usual situation is that the thermal pop- 
ulation of excited states in a parent nucleus sets up 
connections to many states in a daughter by the 
GT operator. However, our current knowledge on 
GT transitions from excited nuclear states is very 
poor, and in many cases, it must rely on theoretical 
calculations. 

• The calculation of forbidden transitions involves 
nuclear transitions between different harmonic os- 
cillator shells and thus requires multi-shell model 
spaces. Such a calculation is not feasible for most 
of the conventional shell models working in one- 
major shell bases. The PSM is a multi-shell shell 
model. This feature is desired particularly when 
forbidden transitions are dominated. 

• Isomeric states belong to a special group of nuclear 
states because of their long half-lives. The exis- 
tence of isomeric states in nuclei could alter signif- 
icantly the elemental abundances produced in nu- 
cleosynthesis. There are cases in which an isomer 
of sufficiently long lifetime can change the paths of 



reactions taking place and lead to a different set of 
elemental abundances [iH, IT^ | . It has been shown 
that the PSM is indeed capable of describing the 
detailed structure of isomeric states ppj] . 

The paper is organized as follows. In Sec. II, we briefly 
introduce the PSM concept and describe how shell model 
diagonalization is carried out in the angular-momcntum- 
projcctcd bases. Technical details of calculation of tran- 
sition matrix elements in the projected bases are given 
in Sec. III. Our first example of a GT transition calcula- 
tion is illustrated in Sec. IV, where we first validate the 
model by comparing the calculated structure properties 
(energy levels and electromagnetic transitions) with ex- 
periment. The obtained GT transition results are then 
discussed and the effects caused by the Gamow- Teller 
residual forces and by configuration mixing are studied. 
Finally, the work is summarized and an outlook on future 
applications is given in Sec. V. 

II. OUTLINE OF THE MODEL 

The Projected Shell Model [13 works with the follow- 
ing scheme. It begins with the deformed Nilsson single 
particle basis, with pairing correlations incorporated into 
the basis by a BCS calculation for the Nilsson states. 
The Nilsson-BCS calculation defines a deformed quasi- 
particle (qp) basis. Then the angular-momentum (and 
if necessary, also particle-number, parity) projection is 
performed on the qp basis to form a shell model space in 
the laboratory frame. Finally a two-body Hamiltonian is 
diagonalized in this projected space. 

Let |$) be the qp vacuum and aj, and a\ the qp cre- 
ation operators, with the index v (it) denoting the neu- 
tron (proton) quantum numbers and running over se- 
lected single-qp states for each configuration. The multi- 
qp configurations of the PSM are given for the four kinds 
of nuclei as follows: 



e-e: {|$} , a^t . |*> , |4>) , .a^ |$) , • • •}, 

odd-:/: {at ( |$) ,4 4 4,4* |$) , aj 4 4 .a\ k |$) , 4,4,4* 4,4 m I*) >■••}, 

odd-^ : {4 t |$) ,4,4,4* I*) .4,4,4* I*) ,44,4*44™ !*>,•■•}, 

o-o: {4,4, 1$) , 4,4,4* 4, l $ > . 44,4*4, I*) -4 i 4 j 4 fc 4,4 m 4„ 1$) >•■ •}■ 



Note that the PSM works with multiple harmonic- 
oscillator shells for both neutrons and protons. The in- 
dices v and 7r in Q are general; for example, a 2-qp 
state can be of positive parity if both quasiparticles i 
and j are from the major iV-shells that differ in N by 
AiV = 0, 2, . . . , or of negative parity if i and j are from 
those TV-shells that differ by AiV = 1,3,.... In bases 
{JU, "• • • " denotes those configurations having a higher 



order of quasiparticles, which may practically be ignored 
in the present study. On the other hand, the bases JTJ 
can be easily enlarged if necessary. If the configurations 
denoted by "■ ■ • " are completely included, one recovers 
the full shell model space written in the representation 
of qp excitation. 

The shell model basis states can then be constructed 
by the projection technique. Without losing generality, 
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the PSM wave-function can be written as 

mM)=J2flK« P MK\®n), 



Kk 



where |$ re ) denotes the qp-basis given in (1), and 
21+ 1 



P 1 — 



8tt 2 



dnDi IK (n)R(n) 



(2) 



(3) 



is the angular momentum projection operator [2l| . In 
Eq. J21 , D MK is the D- function [2^ , R the rotation op- 
erator, and SI the solid angle. If one keeps the axial sym- 
metry in the deformed basis, D MK in Eq. reduces to 
the small d-function and the three dimensions in Q re- 
duce to one. The energies and wave functions (given in 
terms of the coefficients ff K in Eq. J5J) are obtained by 
solving the following eigen-value equation: 



^ ( H Kk,K>k> ~ E I N K K ,K'n')flK K , - 



(4) 



K' 



where H i Kk k , k , and Nj (K k , k , are respectively the matrix 
elements of the Hamiltonian and the norm 



>«> = (®«\Pkk>\®*')- 



(•5) 



The PSM uses a large single-particle space, which en- 
sures that the collective motion and the cross-shell in- 
terplay are defined microscopically by accommodating 
a sufficiently large number of active nucleons. It usu- 
ally includes three (four) major harmonic-oscillator shells 
each for neutrons and protons in a calculation for de- 
formed (superdeformed or superheavy) nuclei. However, 
the shell model dimension in Eq. (0 is small. This means 
that each of the configurations in |J2J is a complex com- 
bination of spherical shell model basis states. Although 
the dimension where the final diagonalization is carried 
out is small, it is huge in terms of original shell model 
configurations. In this sense, the PSM is truly a shell 
model in a truncated multi-major-shell space. 

The Hamiltonian in the present study consists of the 
separable forces 



H = H + H t 



QP 



H, 



GT, 



(6) 



which represent different kinds of characteristic correla- 
tions between valence particles. It has the single-particle 
term Hq 7 the quadrupole+pairing force Hqp with inclu- 
sion of the quadrupole-pairing term, and the Gamow- 
Teller force of the charge-exchange terms Hqt- Ho con- 
tains a set of properly adjusted single-particle energies in 
the Nilsson scheme [23|. The second force, Hqp, contains 



three terms 0| 



Hqp — 



(7) 



which are quadrupole-quadrupole, monopole-pairing, 
and quadrupole-pairing interactions, respectively. The 
strength of the quadrupole-quadrupole force XQQ is de- 
termined in a self-consistent manner that it would give 
the empirical deformation as predicted in mean-field cal- 
culations 0| • The monopole-pairing strength is taken 
to be the form 



Gm — 



G\ =F G2 



N-Z 



A 



(8) 



where "+" ("— ") is for protons (neutrons), N, Z, and A 
are respectively the neutron number, proton number, and 
mass number, and G\ and G2 are the coupling constants 
adjusted to yield the known odd-even mass differences. 
The quadrupole-pairing strength Gq is taken to be about 
20% of Gm, as is often assumed in the PSM calculations 
[l^ . It has been shown in many previous publications 
that such a set of interaction can reasonably well describe 
structures in heavy nuclei. The one-body operators (for 
each kind of nucleons) in Eq. JJJ are of the standard form 

Q2^ = ^2 (a\Q 2fl \a') c* a c a > 

a,a f 
a 

p i = In (aiQ 2 >'>44 

a, a' 

where c£, is the nucleon creation operator, with a stand- 
ing for the quantum numbers of a single-particle state in 
the spherical basis (a = {n£jm}). The time reversal of 
c a is defined as c a = (-) J ~ m c n ^-_ m . 

The last force, Hqt in Eq. ©, is the Gamow- Teller 
force 

H GT = +2 XG T^/3 1 -(-im% 

- 2K GT j2K(- 1 ^ t f-»- ( 9 ) 

This is a charge-dependent separable interaction with 
both particle-hole (ph) and particle-particle (pp) chan- 
nels, which act between protons and neutrons. This t ype 
of force has been used by several authors |24l I25H26I |27| 
in the study of single- and double-/? decay. The ph and 
pp interactions are defined to be repulsive and attractive, 
respectively, when the strength parameters xgt and kqt 
take positive values. The pp interaction, which was intro- 
duced by Kuz'min and Soloviev |24[, is a neutron-proton 
pairing force in the J 77 = 1 + channel. In the present 
calculation, we adopt the interaction strengths 



Xgt = 23/A 
kgt = 7.5/A, 



(10) 



from the original work of Kuz'min and Soloviev [24j. We 
notice that there are different versions for these param- 
eters, for example, the one extracted from systematical 
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studies by Homma et al. |25|. The one-body operators 
appearing in Eq. @ are defined as 

K = T,„,uW\^T-\v)cic v , ^ = (-)"M 

K = f+ = (-)"(rr_ p )t, 

(ii) 

where <r and r are the Pauli spin operator and the isospin 
operator, respectively. As one will see later in an exam- 
ple of the 164 Ho — > 164 Dy transition, the Gamow- Teller 
force is important for a correct reproduction of the GT 
strengths. 

We finally mention that the Hamiltonian assumed in 
the form of Eq. JfjJ may need to be extended when spe- 
cific quantities or transition processes are studied. For 
example, the spin-dipole force would be necessary to re- 
produce first-forbidden transitions. 



where the operators $^ are those given in Eq. A 
GT transition involves a transfer of one unit of angular 
momentum, which is described by the Pauli spin opera- 
tors o^. The isospin operators t± transform a neutron 
into a proton, and vice versa. 

The heart of the present development is the evaluation 
of transition matrix element in the angular-momentum- 
projected bases. Let us start with the PSM wave function 
in Eq. J2J 



Mm)=Y.^P I mk\^ k )- 

Kk 



III. TRANSITION MATRIX ELEMENTS IN 
THE PROJECTED STATES 

The probability of Gamow- Teller transition from an 
initial state of spin Ii to a final state If is defined as 

B{GT,h - I f ) = ^±l\^ If \\^ It )\^ (12) 



In a general weak interaction process with isospin ex- 
change, the initial and final states must correspond to 
different nuclear systems. Therefore, two different sets 
of quasiparticle are generally involved. For example, in 
a calculation of /3-decay from an odd-odd parent to an 
even-even daughter nucleus, |$ K ) for the initial odd-odd 
system has the form 



|$«(a)) = {at a* \a) , ata£4 a t |a) , 4444 \a) , 444444 \a) ,■■■}, (13) 
and for the final even-even system, 

!$„(&)) = {\b) , 6t ft t | 6) ; 6 t 6 t \ b) 1 6 t 6 t 6 t 6 t | 6) 7 . . . } . (14) 

Here, {a^} and {b^} are two different sets of quasiparticle operators associated with the quasiparticle vacua \a) and 
\b), respectively. 

To compute the matrix element of a tensor operator of rank A in the projected states, the matrix element of the 
operator Pll M Tx^P^^, has to be evaluated. A useful expression was derived in Ref. [l2| 



K/M/ 1 ^ r MiKi 



(15) 



which follows from the transformation property of a tensor operator of rank A under rotation as well as the reduction 
theorem of a product of two D-functions [23] ■ With the help of this relation, the reduced matrix element in Eq. i|12|) 
can be written as 



KiK,KfK,' (. n ) 



(16) 



If we assume an axial symmetry in nuclei, as for the 
example we shall discuss in the next section, the gen- 
eral three-dimensional angular momentum projection is 
reduced to a problem of one-dimensional projection, with 



the projector having the following form 



PLk = I /+ - 



d/3 sin/3 <f MK (fi) R y (p) (17) 
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with 



Ry(j3) = e" jpj ". (18) 
In Eq. HI 7(1 > <*mk (/?) is the small-d function and j3 is one 



J 



of the Euler angels Inserting the projector 1 (17(1 into 
(^'WI/^^Kz-ntfJ^C )) in E q- ©> one obtains 



(*«<(b)\P±P I J Ki \$ K (a)) = 



d/3 sin/3 ^W^K^'WI/tA^I^M) 
I 



(19) 



The problem is now reduced to evaluation of the matrix 
element 



($ K ,(b)\/3tMf3)\* K (a)}, 



(20) 



which is the problem of calculating the (5 operator sand- 
wiched by amulti-qp state | $«/(&)) and a rotated multi-qp 
state R y (/3) |$ K (a)), with a and 6 characterizing different 
qp sets. 

To calculate /3 X Ry{0) |$ K (a)), one must com- 

pute the following types of contractions for the Fermion 
operators 



Ai 
B 



'13 n-Hfj-j 

where we have defined 



(b\ \p]a\a]\a) = [V{(3)U- 1 (P)h, 
(b\bMa]W) = [U- 1 (P)] l] , 



(21) 



(6|i? y (/3)|«) 



and 



(b\R y (/3)\a) = [det f/(/3)] 1/2 . 



(22) 



Eqs. ((21(1 and 1(22(1 are written in a compact form of 
N x N matrix, with N being the number of total single 



particles. The general principle of finding U(f3) and V((3) 
is given by the Thouless theorem [2^ , and a well worked- 
out scheme can be found in the work of Tanabe et al. 

m. 

To write the matrices U(/3) and V{f3) explicitly, we 
consider the fact that {ai,a\} and {bi,b\} can both be 
expressed by the spherical representation {q, cj} through 
the HFB transformation 



c 




V a 


ct_ 


- u« 


u a 


c 




v b 


c t_ 


^ v b 





(23) 



U a ,V a ,Ub and Vf, in above equations, which define the 
HFB transformation, are obtained from the Nilsson-BCS 
calculation. A rotation of the spherical basis can be writ- 
ten in a matrix form as 



Ry(0) 







(24) 



Combining Eqs. ((23(1 and ((24(1 and noting the unitarity 
of the HFB transformation, one obtains 



Ry{(3) 



b 



d[fi) 


1 


( Ua 


M 


a 









U a J 





C/(/3) and can finally be obtained from the following equation 

( d((3) 



(U(J3) V({3) 
{ V(P) U(f3) 



v b T ul )\ d(P) 



U a V a 
V a U a 



U^d((3)Ua + V b T d(f3)V a U^d((3)V a + V? 'd{(3)U a 
V b T d{(i)U a + U^d{(3)V a U^d{f3)U a + V? d(0)V a 
I 



(25) 



(26) 



We identify the following symmetry property for the 



matrix element 1(20(1 



\ "* jl 



dl.-M(*«\Pw R (0)\*«')- 



(27) 
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In addition, there is a symmetry relation for the reduced 
matrix element 



(*Z / ||)9 ± ||*7 < > = H 



2L + 1 



21,+ 



I (^ i ||/3 = F||*z / ).(28) 



These symmetry properties can be used for testing the 
coding and checking the numerical accuracy because the 
matrix elements at each side of the equations are calcu- 
lated independently. 



IV. GAMOW-TELLER TRANSITION RATES IN 
A RARE-EARTH EXAMPLE 

The formalism derived in the preceding section is valid 
for general /3-decay and electron-capture calculations no 
matter whether a process is associated with an allowed or 
a forbidden transition. Our following illustrative example 
considers allowed transition only, for which the following 
selection rules (change in quantum numbers between the 
initial and the final state) apply 

AI = 0,±1, 
An = +1. 

We take the odd-odd nucleus 164 Hog7 as the example, 
which decays to 164 Dy 98 through an electron-capture pro- 
cess by changing a proton to a neutron. We shall evaluate 
B(GT) values defined in Eq. (TgJ. For the A = 164 nu- 
clei, the GT interaction strengths in Eq. (|10|l become 

Xgt = 23 /A = 0.14 (MeV) 
KGT = 7.5/A — 0.046 (MeV). 

In order to compensate for the dependence of the decay 
rate on the transition energy, it is customary to express 
the transition probability in terms of the product ft |30j| , 



6163.4 



te) e ff^ GT ) 



(29) 



where t is the half-life and / is a dimensionless quantity 
depending on the charge of the nucleus and the energy 
and multipolarity of the transition. For the coupling con- 
stants appearing in Eq. (}29§ , we simply adopt those from 
Ref. |3j| (see also the references cited therein), 



9 A =0.74^, 
9vJeS 9v 



M = -1.26. 

gv 



A. Structure description: Energy levels and 
electromagnetic transitions 

We use the standard Nilsson scheme [23| to generate 
deformed single particle states for our basis. For the Nils- 
son parameters we use those of Jain et al. |32| , which are 
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FIG. 1: (Color online) Calculated energy levels of the ground 
band in 164 Ho (upper panel) and 164 Dy (lower panel) are com- 
pared with data. The ground state of the odd-odd 164 Ho 
nucleus has I n — 1 + . 



1 S f 



3 


400- 








300- 


1 




1 


200- 


cm" 




LU 


100- 


CD 





164 



Dy 



• Exp. 
-□-PSM 



H — • — r - 

2 4 



10 12 



Spin (h) 



FIG. 2: (Color online) Calculated B(E2) values of the ground 
band in 164 Dy are compared with data. 



applicable to the rare earth region. The deformed bases 
are obtained with the quadrupole deformation £2 = 0.28. 
At this deformation, the Fermi levels of the A = 164 nu- 
clei are surrounded by several neutron and proton single 
particle orbitals, such as v [523] f and tt[523]|, which are 
the most relevant ones in our later discussion on GT- 
transition. The monopole-pairing coupling constants in 
Eq. © are taken to be G x = 20.12 MeV and G 2 = 13.13 
MeV, which are the same values employed in many pre- 
vious rare earth calculations. The quadrupole-pairing 
strength Gq is proportional to Gm, the proportionality 
constant being fixed to be 0.20 for the nuclei considered 
in this paper. 

Under the above calculation conditions, energy levels 
and wave functions are obtained for the parent 164 Ho 
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and the daughter 164 Dy. In Fig. 1, we compare the 
calculated energy levels of the ground bands with data. 
An excellent agreement between theory and experiment 
is achieved, implying that with this established set of 
parameters, the structure of these nuclei can be micro- 
scopically described. In Fig. 2, one sees that the calcu- 
lated E2 transition probabilities along the ground band 
in 164 Dy also agree well with data. These B(E2) values 
are calculated with the resulting wave functions as 



B{E2Ji^I f ) 



2J/ + 1 
21, + 1 



(* If \\Q 2 \\y u ) 



(30) 



where the effective charges 1.5e for protons and 0.5e for 
neutrons are used. 
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FIG. 3: Distribution of the B(GT+) strength. 



The above calculations as well as those of the previ- 
ous PSM publications indicate that for low-lying states 
in heavy, deformed nuclei, where one has a lot of struc- 
ture data to compare with, the PSM can usually pro- 
vide a good description. For GT transitions, a com- 
plete description requires the model to cover those highly 
excited states as well because the GT strength usually 
spreads for a wide range of excitation and a large amount 
of GT strength can concentrate in some high energy 
states. Next, we examine the resulting wave functions 
of both low-lying and high-lying states in the full model 
space. There is a well-known charge-exchange sum-rule, 
the Ikeda sum-rule 33] , expressed as 

S(GT-) -S'(GT + ) 

= 5>(GT-,z^/)-^5(GT+,z^/) 
/ / 



f,li 

3(N-Z). 



/■/' 



(31) 



This model independent sum-rule results from the com- 
mutation relation between the isospin operators. Hence, 
a (non-relativistic) model is expected to satisfy it as long 
as the basis in Eq. (|3*TJl is complete. 

Let in Eq. (|3l")l correspond to the 0-qp ground 
state of a nucleus and \^ff) to the f-th 2-qp state with 
I* = 1 + . The summation in Eq. (|31ll runs over all pos- 
sible 2-qp 1 + states in the complete model space. When 
the single-particle states with N — 3, 4, 5, 6 harmonic- 
oscillator shells are taken into account for both neutrons 
and protons, the total number of 1 + states turns out to 
be 1718. For our example of 164 Dy (N = 98 and Z = 66), 
we show in Figs. 3, 4, and 5 all the obtained B(GT + ) and 
B(GT~) strengths as functions of the 1 + state energy. In 
Fig. 3, it is observed that a considerable portion of the 
B(GT + ) strength appears at the low excitation region 
less than 2.5 MeV. The rest of the strength is fragmented 
over many states extended to about 20 MeV. The frag- 
mentation is clearly caused by the residual interactions. 
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FIG. 4: Distribution of the B(GT~) strength. 
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FIG. 5: Distribution of the B(GT~) strength. The numbers 
are the same as those of Fig. 4, but plotted in the scale as 
Fig. 3. 
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In Fig. 4, a strong peak of the B(GT~) strength appears 
around 16 MeV, while the rest of the strength spreads 
for nearly the entire energy region. The concentration 
around 16 MeV takes most of the B(GT~) strength, and 
has a character of GT giant resonance. In Fig. 5, we show 
the same B(GT~) distribution as that of Fig. 4, but with 
a different scale to emphasize those smaller strengths. It 
is seen that the fragmentation extends to the excitation 
of 30 MeV. 

Adding respectively all the B(GT+) and B(GT~) 
strengths together, as indicated in Eq. we get 

S(GT-) = 97.4957, S(GT + ) = 4.1717. 

This means that our 5(GT~) - S(GT+) exhausts 97.2% 
of the sum-rule value 3(N — Z). Note that in our Nilsson 
calculation for deformed single particle states, we use the 
standard formalism 23] in which the neutron harmonic- 
oscillator frequency lo v differs from the proton frequency 
uj n for N ^ Z nuclei. If we set ui v = u w , the calculated 
values become 

S(GT-) = 97.7118, S(GT+) = 1.7797, 

which approaches 99.9% of the sum-rule value. If we fur- 
ther switch off all the pairing interactions, i.e. if we aban- 
don the quasiparticle picture and go back to the particle 
picture, we get 



S(GT-) = 96.7480, S(GT^ 



0.7480. 



Thus, in the idealized (but not realistic) case in which 
neutrons and protons move in the same harmonic- 
oscillator potential well with no pairing interaction, the 
Ikeda sum-rule 1)31(1 is fulfilled exactly. The above results 
may be taken as an additional nontrivial check for the 
coding since in all these calculations, the wave functions 
in Eq. have a form of Eq. © and the B(GT) val- 
ues are computed in a sophisticated way as presented in 
Sections II and III. 



C. Effects of the Gamow- Teller forces and 
configuration mixing 

The GT interaction in Eq. is a new addition to 
the PSM Hamiltonian ©. The results in Figs. 1 and 2 
seem to indicate that the GT forces do not have effects 
on the low-lying states since with or without the GT 
forces, the energy levels and electromagnetic transitions 
are described equally well. The purpose of this section is 
to study how the new addition of the GT forces influences 
the B(GT) values and in what extent the configuration 
mixing modifies the results. 

We perform calculations for the GT transition from 
the 1 + ground state of 164 Ho to the + ground state of 
164 Dy with five different cases as listed in Table I, and 
results are shown in Fig. 6. For Cases 1 and 2, we do not 
include the GT forces. In Case 1, there is only one di- 
mension in both the parent 164 Ho and the daughter 164 Dy 



TABLE I: Five different cases of the GT transition calcula- 
tion. In the table, "even-even" means the + ground state of 
164 Dy and "odd-odd" means the 1 + ground state of 164 Ho. 
Hgt stands for the Gamow- Teller forces in the Hamiltonian. 

Case I Even-even I Odd-odd I Hgt 



{|0» (1) 

{|0>, \2v2-k)} (857) 

{|0>} (1) 

{JO), \2u2-k)} (322) 

{|0>, |2^2tt}} (857) 



•H523] 



x^[523]|} (1) 

{|1zTtt)} (823) 

{jlflTr)} (823) 

{|1j/1tt)} (823) 

{jli/lvr)} (823) 



no 
no 

yes 
yes 
yes 



0.8 



0.6 



CD 0.4 

CO 



0.2- 



0.0 




H 1 

2 3 
Case 



FIG. 6: (Color online) B(GT) values calculated for the five 
different cases. Available data are shown in stars. The B(GT) 
data are converted from the experimental Log/t values 
through Eq. 1^51 . 



wave function; namely, we approach the wave function of 
the even-even nucleus using the projected vacuum state 
and that of the odd-odd nucleus using the projected 2- 
qp state {z/[523]§ x tt[523]|}. In Case 2, we include 823 
2-qp basis states in the odd-odd and additional 856 4-qp 
states in the even-even wave function. As one can see 
from Fig. 6, there is no much difference between Case 1 
and Case 2, implying that without the GT forces, con- 
figuration mixing becomes unimportant. It is also seen 
in these two cases that without the GT forces, the calcu- 
lated B(GT) is by nearly a factor of three too large when 
compared with data. 

A significant reduction in B(GT) is obtained when we 
switch on the GT forces. As shown in Fig. 6 for Cases 3, 
4, and 5, these forces push down the B(GT) values con- 
siderably. Moreover, effect of configuration mixing shows 
up when the GT forces are included; there is a clear dif- 
ference between the results of Case 3 (with no 4-qp state 
in the even-even wave function) and Case 4 (with 321 
4-qp basis states mixed in the even-even wave function) . 
The latter case reproduces data well. Comparing the re- 
sult of Case 5 with Case 4, one sees a converging trend in 
the calculation; further enlarging the configuration space 
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TABLE II: Comparison of calculated Logft values with data 
for electron-capture from 164 Ho to 164 Dy and /3~ decay from 
164 Ho to 164 Er. Data are taken from Ref. H. 



Decay 



iLog/i (Th.) Log/t (Exp. 



164 Ho (1+) 
164 Ho (1+) 



164 Dy (0+) 
164 Dy (2+) 



4.61 
4.96 



4.59 (8) 
4.86 (6) 



164 Ho (l 4 
164 Ho (V 



164 Er (0+) 
164 Er (2+) 



5.19 
5.62 



5.54 (10) 
5.75 (10) 



1.6- 
1.4- 
1.2- 
> 1.0- 
^;0.8- 
w 0.6- 
0.4- 
0.2- 
0.0- 



8 T 

— Exp. r P 

-□-PSM K *P 

-□ /164 H0 




"i — 1 — r 



2 4 
Spin (Pi) 



FIG. 7: (Color online) Possible transition paths from the 
ground band of 164 Ho to that of 164 Dy. Experimentally known 
energies for these states are also displayed. 



would no longer improve the result as far as the GT tran- 
sition between the ground states is concerned. 

164 Ho can capture an electron and change to 164 Dy, 
and can /3~-decay to 164 Er. In Table II, we compare 
calculated Logft values with the existing data for both 
processes. We have used the same parameters in gen- 
erating the 164 Er wave functions and in the calculation 
of the transition matrix elements. As can be seen from 
Table II, the agreement in both processes is satisfactory. 



D. GT transition rates of the excited states 

In Fig. 7, possible allowed GT transitions from the 
ground band of 164 Ho to that of 164 Dy are illustrated. 
Here, energy levels up to about 800 keV of excitation in 
both parent and daughter nuclei are considered. They are 
AI = I p — I c i = +1 transitions (in black), AI = transi- 
tions (in red), and AI = — 1 transitions (in green). The 
calculated B(GT) values are shown in the upper panel, 
and the corresponding Logft values in the lower panel of 



0.3 -. 

^ 0.2 - 
h- 

0.1 -| 

0.0 



CQ 




"I 1 1 1 1 1 1 1 

2 4 6 8 

Spin I (h) 
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4.5 - 
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• 




i 
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Spin I (fi) 
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FIG. 8: (Color online) Calculated B(GT) (upper panel) and 
Logft values (lower panel) for the 164 Ho — > 164 Dy electron- 
capture process. Available data are shown in stars. Data are 
taken from Ref. 1341 . 



Fig. 8. The experimentally measured decay probabilities 
are those of the I p = 1 — » Id = and I p = 1 — » Id = 2 
transitions, with which our calculation agrees well. The 
other transition probabilities associated with the excited 
states are our prediction. Note that the decay rates with 
AI = (in red) are predicted to have larger B(GT) and 
smaller Logft values than the measured AI = ±1 transi- 
tions. This is an interesting result because normally, it is 
very difficult to study /3-decay rates of excited states in 
laboratory since excited nuclear states decay much faster 
via 7-emission than by /3-decay. A measurement may 
be possible only if the underlying states are long-lived 
isomers. 

In stellar environments, nuclei that are exposed to high 
temperatures and high densities may experience a signif- 
icant enhancement in their decay rate. This idea was 
brought out early by Cameron [35| and studied in detail 
by Bahcall 36] . The enhancement can result from the ef- 
fects related to the thermal population of excited nuclear 
states in the hot photon bath. In a thermal equilibrium 
with temperature T, the population probability is deter- 
mined by the Boltzmann factor and the statistical weight 



Pi 



(2h + l)e- g -/ fcT 
£ m (2/ m + l) e -JW*T- 



(32) 



Hence the actual stellar decay rate must consider decays 
of all thermally populated excited states i in the parent 
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nucleus to the accessible levels j of the daughter. Ther- 
mal enhancement of beta-decay rates can have substan- 
tial impact on nucleosynthesis, for example, the s-process 
branching conditions [37j ■ 

V. SUMMARY AND OUTLOOK 

In this article, we have presented the new development 
of a shell model method for calculation of Gamow- Teller 
transition rates. The method is based on the Projected 
Shell Model. Different from the conventional shell model, 
which builds its states in a spherical basis, the PSM con- 
structs its states in a deformed basis in which impor- 
tant nuclear correlations are taken into account very effi- 
ciently Therefore, it is possible for a shell model diago- 
nalization in the PSM to be carried out in a manageable 
space for medium to heavy, and even for super-heavy nu- 
clei. 

We have shown how the GT transition matrix elements 
are calculated in the PSM framework. A computer code 
has been developed and been tested. One nontrivial test 
has been done through the Ikeda sum-rule. We have ob- 
tained a reasonable distribution of the B(GT) strength 
with a fulfillment of the sum-rule. We have presented the 
first example from the rare earth region. In the calcula- 
tion of 164 Ho — > 164 Dy electron-capture process, we have 
predicted the GT transition rates for the excited states. 
Such rates should be included as part of the total rate 
when these states are thermally populated in hot stellar 
environments. 

In the calculation discussed so far, we have considered 
only allowed /3-decay for the low-lying states. Calculation 



of allowed /3-decay for states with high excitation and the 
forbidden transitions is possible in the PSM framework. 
Study of GT giant resonance is also under consideration. 

The method described in the present article can be 
applied to various fields such as nuclear astrophysics and 
fundamental physics, where weak interaction processes 
take place in nuclear systems j33|. In particular, one may 
find interesting applications to cases where a laboratory 
measurement for certain weak interaction rates is diffi- 
cult and where the conventional shell model calculations 
are not feasible. Potential applications in nuclear astro- 
physics are calculations of /3-decay rates for the r-process 
[39| and the rp-process |4£j nucleosynthesis, and electron- 
capture rates for the core collapse supernova modelling 
0, |4l| . In the double-/? decay theory, theoretical calcu- 
lations for the nuclear matrix elements are needed, for 
which one has relied on the Quasiparticle Random Phase 
Approximation |26t l42j| , particularly when heavy nuclei 
are involved. We expect that the method presented here 
can make important contributions to all these studies. 
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